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ABSTRACT 

We derive constraints on cosmological parameters using the power spectrum of galaxy 
clustering measured from the final two-degree field galaxy redshift survey (2dF- 
GRS) and a compilation of measurements of the temperature power spectrum and 
temperature-polarization cross-correlation of the cosmic microwave background radi- 
ation. We analyse a range of parameter sets and priors, allowing for massive neutrinos, 
curvature, tensors and general dark energy models. In all cases, the combination of 
datasets tightens the constraints, with the most dramatic improvements found for the 
density of dark matter and the energy-density of dark energy. If we assume a flat 
universe, we find a matter density parameter of = 0.237 ±0.020, a baryon density 
parameter of Ob = 0.041 ± 0.002, a Hubble constant of i?o = 74 ± 2 kms"^Mpc"\ a 
linear theory matter fluctuation amplitude of CTg = 0.77 ± 0.05 and a scalar spectral 
index of rig = 0.954 ± 0.023 (all errors show the 68% interval). Our estimate of rig 
is only marginally consistent with the scale invariant value rig = 1; this spectrum is 
formally excluded at the 95% confidence level. However, the detection of a tilt in the 
spectrum is sensitive to the choice of parameter space. If we allow the equation of 
state of the dark energy to float, we flnd wde = — 0.85tQ;^y, consistent with a cosmo- 
logical constant. We also place new limits on the mass fraction of massive neutrinos: 
fi, < 0.105 at the 95% level, corresponding to X)"^'' < 1-2 eV. 

Key words: large scale structure of the universe - cosmic microwave background, 
cosmological parameters 



1 INTRODUCTION 

Since the turn of the millennium, we have witnessed a dra- 
matic improvement in the resolution and accuracy of mea- 
surements of fluctuations in the temperature of the cosmic 
microwave background radiation (CMB). The discovery of 
features in the power spectrum of the CMB temperature, 
the acoustic peaks, marked the start of a new, data-rich era 
in cosmology (de Bernardis et al. 2000; Hanany et al. 2000). 
The relative positions and heights of the acoustic peaks en- 
code information about the values of the fundamental cos- 
mological parameters, such as the curvature of the universe 
or the physical density in cold dark matter and baryons. Per- 
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haps the most striking example of the progress achieved is 
the first year data from the WMAP satellite (Bennett et al. 
2003; Hinshaw et al. 2003). 

The CMB data alone, however, do not constrain all of 
the fundamental cosmological parameters to high precision. 
Degeneracies exist between certain combinations of parame- 
ters which lead to indistinguishable temperature fluctuation 
spectra (Efstathiou & Bond 1999). Some of these degen- 
eracies can be broken by comparing theoretical models to 
a combination of the CMB data and other datasets, such 
as the power spectrum of galaxy clustering. At the same 
time as the new measurements of the CMB were obtained, 
two groundbreaking surveys of galaxies in the local Universe 
were being conducted. The two-degree field galaxy redshift 
survey (2dFGRS; CoUess et al. 2001; 2003) and the Sloan 
Digital Sky Survey (SDSS; York et al. 2000; Abazajian et al. 
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2005) are substantially larger than previous redshift surveys 
and allow the clustering of galaxies to be measured accu- 
rately on all scales. On large scales, the connection to theo- 
retical models is most straightforward. 

Percival et al. (2001) used the power spectrum of galaxy 
clustering measured from the 2dFGRS to constrain the ra- 
tio of the baryon to matter density, nb/f2m, and the matter 
density, Q,ynh. Efstathiou et al. (2002) used a compilation of 
pre-WMAP CMB data and the Percival et al. measurement 
of the galaxy power spectrum to find conclusive evidence for 
a non-zero cosmological constant, independent of the Hub- 
ble diagram of distant Type la supernovae. Percival et al. 
(2002) again used pre-WMAP CMB data and the early 2dF- 
GRS power spectrum measurement to place constraints on 
cosmological parameters in flat models. The WMAP team 
also used the Percival et al. galaxy power spectrum in their 
estimation of cosmological parameters (Spergel et al. 2003). 
Other papers have also analyzed the information encoded 
in the 2dFGRS and SDSS power spectra (Tegmark, Zaldar- 
riaga & Hamilton 2001; Pope et al. 2004; Tegmark et al. 
2004b; Seljak et al. 2005). In view of the impact of this 
work, the recent completion by Cole et al. (2005) of the 
power-spectrum analysis of the final 2dFGRS dataset is an 
important development. The Cole et al. results are nearly 
twice as accurate as those obtained from the part-complete 
2dFGRS in 2001, and a key aim of the current paper is to 
see how this affects the outcome of joint analyses including 
CMB data. 

In view of these rapid improvements in our knowledge 
of the cosmological parameters, it is also important to take 
stock of precisely which parts of the model are actually be- 
ing tested. Quite often, restrictive assumptions have been 
adopted for the background cosmology when claims are 
made about the constraints on a particular parameter. It 
is important to establish how robust the constraints really 
are when the data are compared with more general cosmo- 
logical models. 

Our goal is thus to establish firmly how well the lat- 
est CMB and large-scale structure (LSS) data determine a 
broad set of cosmological parameters, paying attention to 
how the choice of priors for parameter values and the com- 
bination of different parameters can influence the results. 
The outline of the paper is as follows. In Section 2, we de- 
scribe the data used in our parameter estimation and set 
out the various parameter spaces studied. In Section 3, we 
present our main results for the parameter constraints ob- 
tained by comparing theoretical models to the CMB data 
and the galaxy power spectrum of the final 2dFGRS mea- 
sured by Cole et al. (2005). We assess the impact of dif- 
ferent choices for priors and parameter sets in Section 4. 
We explore the justification for using models with different 
numbers of free parameters in Section 5. In Section 6, we ex- 
amine how the parameter constraints change when the SDSS 
galaxy power spectrum measured by Tegmark et al. (2004b) 
is used instead of the 2dFGRS power spectrum. Finally, we 
summarize our conclusions in Section 7. 



2 THE METHOD 

We now set out the approach we will take to constrain the 
values of the basic cosmological parameters. In Section [2. II 
we list the CMB and LSS datasets that we compare against 
the theoretical models and explain how these datasets are 
modelled. The parameter sets that we will consider are de- 
fined in Section f2.2l The methodology for searching param- 
eter space and placing constraints on parameters is set out 
in Section [2.31 

2.1 The datasets 

In order to constrain the parameters in our cosmological 
model, we use a compilation of recent measurements of the 
CMB and the power spectrum of galaxy clustering in the 
local Universe: 

(i) The WMAP flrst year temperature power spectrum for 
spherical harmonics 2 <l < 900 (Hinshaw et al. 2003). 

(ii) Observations of the temperature spectrum over the 
spherical harmonic range 900 < £ < 1800 made up to July 
2002 using the Arcminute Cosmology Bolometer Array Re- 
ceiver (ACBAR; Kuo et al. 2004). 

(iii) The temperature spectrum for 600 < £ < 1500 mea- 
sured using the Very Small Array (VSA; Dickinson et al. 
2004). 

(iv) Two years of temperature correlation data with 600 < 
£ < 1600 from the Cosmic Background Imager (CBI; Read- 
head et al. 2004). 

(v) The WMAP first year temperature-polarization power 
spectrum for spherical harmonics 2 < £ < 450 (Kogut et al. 
2003). 

(vi) The power spectrum of galaxy clustering measured from 
the final 2dFGRS catalogue (Cole et al. 2005). 

The four measurements (i-iv) of the power spectrum of 
temperature fluctuations in the CMB extend over the spher- 
ical harmonic range 2 < £ < 1800. Some of the available 
datasets extend to higher multipoles. However, we do not 
include these scales in our analysis, as the temperature fluc- 
tuations on such scales can be strongly affected by secondary 
sources. The WMAP team adopted a similar approach, aug- 
menting the first year WMAP data with other experiments 
which have better angular resolution (Spergel et al. 2003). 
However, the VSA data were not available to the WMAP 
team at the time that the paper by Spergel et al. was writ- 
ten. Theoretical temperature-temperature and temperature- 
polarization spectra are computed for each model using CAMB 
(Lewis, Challinor & Lasenby 2000). 

Cole et al. (2005) measured the power spectrum of 
galaxy clustering from the final 2dFGRS catalogue. The 
power spectrum measured for galaxies differs in a number 
of ways from the power spectrum for the mass predicted in 
linear perturbation theory; (i) Nonlinear evolution of den- 
sity perturbations leads to coupling between Fourier modes, 
changing the shape of the power spectrum, (ii) The galaxy 
power spectrum is distorted by the gravitationally induced 
peculiar motions of galaxies when a redshift is used to in- 
fer the distance to each galaxy, (iii) The power spectrum of 
the galaxies could be a modified version of the power spec- 
trum of the mass. This phenomenon is known as galaxy bias. 
The ratio between the galaxy and matter spectra could also 
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change with scale. (However, we assume a constant bias over 
the scales considered in this paper), (iv) The power spectrum 
measured by Cole et al. from the 2dFGRS is the direct trans- 
form of the data, and is thus what CMB researchers would 
term a pseudo-spectrum. As such, it yields a convolution 
of the underlying galaxy power spectrum with the modulus 
squared of the Fourier transform of the window function of 
the survey. 

In order to constrain cosmological parameters, these ef- 
fects need to be modelled. The accuracy of the modelling 
requires that the comparison between theory and observa- 
tion should be restricted to a limited range of scales. We use 
the 2dFGRS power spectrum data for k < 0.15 /iMpc~^ and 
discard measurements with k < 0.02 /iMpc~^ which could 
be affected by uncertainties in the mean density of galax- 
ies within the survey. We follow the scheme used by Cole 
et al. who applied a correction for non-linearity and scale- 
dependent bias to the shape of P{k) of the form 

P.Ak)=b'^^^Pun{k), (1) 

where A = lA and Q — 4.6 are the preferred values and b is 
a constant bias factor. This formula is deduced by compari- 
son with detailed numerical galaxy-formation models: these 
show that the value of A is robust, but the exact value of 
Q depends on galaxy type and also has some uncertainty 
depending on how the modelling is done. These results were 
used to determine a range of allowed values for Q, from 
which the value Q — 4.6 is preferred; with this choice, robust 
parameter constraints are obtained if one considers maxi- 
mum k values beyond our limit of 0.15/iMpc~^. For this 
limit, neglecting the correction entirely and simply fitting 
linear theory yields almost identical results to those pre- 
sented here. In particular, it has no impact on the marginal 
indication of a deviation from ris = 1. 

2.2 The parameter space 

In this paper, we make the basic assumption that the pri- 
mordial density fluctuations were adiabatic, Gaussian and 
had a power law spectrum of Fourier amplitudes. As pointed 
out by Leach & Liddle (2003a), the CMB data prior to 
WMAP were of insufficient quality to justify the rejection of 
this simple hypothesis. Following the release of the WMAP 
first year results, which do have the precision required to 
test this model, our assumptions remain well motivated. Ko- 
matsu et al. (2003) found that the WMAP sky maps are 
consistent with Gaussian primordial fluctuations to a much 
higher precision than was attainable with COBE. Peiris et 
al. (2003) found that models with a spectral index vary- 
ing slowly with wavenumber give slightly better fits to the 
WMAP data, particularly when combined with estimates of 
the power spectrum of the Lyman-a forest. However, the 
evidence for a running spectral index is weak and has been 
disputed by other groups (e.g. Bridle et al. 2003b; Slosar, 
Seljak & Makarov 2003; Seljak et al. 2005). Bennett et al. 
(2003) and Spergel et al. (2003) point out that, on large 
scales, a few modes of the CMB temperature power spec- 
trum measured by WMAP lie below the predictions of the 
standard ACDM model. One interpretation of this apparent 
discrepancy is that new physics may be needed (e.g. Bridle 
et al. 2003a; Efstathiou 2003). However, several studies have 



argued that the disagreement is actually less significant than 
was first claimed (Gaztanaga et al. 2003; de Oliveira-Costa 
et al. 2003; Efstathiou 2004). 

^From the above starting point, the cosmological model 
we consider is defined by eleven parameters: 

P = (r2fe,(.Jdm,t^b,,f>^,ifDE,T, ris, As,r, 6, e). (2) 

There are eight further basic quantities whose values can be 
derived given the above set: 

= {QnE,h,Qm,crs,Zrc,to,'^2m„,nt). (3) 

We now go through the parameters in these lists, defining 
each one and explaining how the values of the derived pa- 
rameters are obtained. 

The are five quantities that describe the homogeneous 
background cosmology through various contributions to the 
mass-energy density. These are, in units of the critical den- 
sity: fifc, which describes the curvature of the universe; floE, 
the energy-density of the dark energy; ujdm = ildmh^ , the 
density of the dark matter (where f2dm = ^cdm + i^u is the 
sum of the cold and hot dark matter components and h is 
Hubble's constant in units of 100 kms~^Mpc~^); Ub = ^bh^, 
the baryon density; and fv = i^u/ildm, the fraction of the 
dark matter in the form of massive neutrinos. The sum of 
neutrino masses is given by ^m^. The matter density pa- 
rameter is given by Q^n ~ f^dm -I - fib ■ The value o f the Hubble 
constant is derived from h = a/ (oidm + u^b ) /^m ■ The energy 
density of the dark energy is set by Que = 1 — — ^k- 
The dark energy component is assumed to have an equation 
of state that is independent of redshift, with the ratio of 
pressure to density given by wue- 

There are four quantities that describe the form of the 
initial fiuctuations; the spectral indices, Us and rit, and the 
primordial amplitudes. As and rAs, of scalar and tensor fiuc- 
tuations respectively. These parameter values are quoted at 
the 'pivot' scale wavenumber of fc = 0.05 Mpc~^. We can 
translate the results obtained for As into a constraint on 
the more familiar parameter as, the rms linear perturba- 
tion theory variance in spheres of radius 8 /i~^Mpc, using 
the matter fiuctuation transfer function. Note that when we 
consider tensor modes, we make the slow-roll assumption 
that Tit = —r/S. 

The bias factor, b = ^ Pgai(fc)/-PDM(fc), describes the 
difference in amplitude between the galaxy power spectrum 
and that of the underlying dark matter. The value of 6 is 
marginalized over, using the analytic expression given in 
Appendix F of Lewis & Bridle (2002). We assume that the 
reionization of the neutral intergalactic medium occurred in- 
stantaneously, with an optical depth given by r; the redshift 
of reionization, Zrc, depends upon a combination of param- 
eters (see Table 1 of Tegmark et al. 2004b). The age of the 
universe is to. 

Finally, O gives the ratio of the sound horizon scale at 
the epoch of decoupling to the angular diameter distance to 
the corresponding redshift and replaces the Hubble constant 
as a base parameter (Kosowsky, Milosavljevic & Jimenez 
2002). We have chosen to use this parameter, rather than, 
for example, the energy density in dark energy, since it has a 
posterior distribution that is close to Gaussian. This reduces 
degeneracies between parameters and results in a faster con- 
vergence of our search of parameter space (see Section [2.311 . 
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compared with studies in which parameters such as SIde, 
which does not have a Gaussian posterior distribution, are 
allowed to vary. This approach is a standard practice even 
though, usually, the final results are expressed in terms of 
more familiar parameters such as JIde or h. However, care 
must be taken when comparing our results with those from 
studies which have assumed flat priors on different parame- 
ters in their Bayesian analysis. Such choices may affect the 
final results in a way that is difficult to quantify. 

We do not attempt to vary all eleven parameters of the 
model at once. Such an approach would lead to a mixture 
of poor estimates of the values of individual parameters and 
constraints on various combinations of parameters. As we 
are primarily interested in deriving the best possible con- 
straints on individual parameters, we instead consider sub- 
sets of parameter space, varying five, six or seven parame- 
ters at a time. Of the remaining parameters, some are held 
at fixed values and the others are referred to as derived pa- 
rameters. The values of the derived parameters follow from 
the values of other parameters, once the assumptions made 
in each parameter space have been taken into account. We 
will now set out each of our parameter spaces in turn, stat- 
ing which parameters are varied and which are held fixed. 
In all cases, the bias parameter, h, is marginalized over, so 
we do not include this in the list of parameters whose values 
are constrained. 

In the simplest case, we vary five parameters which we 
refer to as the 'basic-five' (b5) parameter set. The following 
parameters are allowed to float: 



(tJdm,C^b,T, As, 9). 



(4) 



The values of the fixed parameters in the basic-five model 
are: 



Pflxcd = (^^k = = Oj-u^DE = -l,ns = l,r = 0). 



(5) 



The results of this model are discussed in Section 3.2. 

The basic-five set is expanded to allow the value of the 
scalar spectral index to float, giving the basic-six (b6) model 
(see Section 3.3): 

Pvaricd = (t^dm , (^b , T, Tls , As , 9) . (6) 

The fixed parameters in the basic-six model are: 

PL,d = (!^k = 0, = 0, WDE = -1, r = 0). (7) 

We also consider four parameter spaces in which one 
additional parameter is constrained along with the basic- 
six set. In Section 3.3, the additional parameter is the mass 
fraction of massive neutrinos, /^: 

Pvifcd = (t^dm, <^b, fu, T, Tls, As, 9). (8) 



The fixed parameters in this case are: 
Pfiid = («k = 0,™de = -l,r = 0). 



(9) 



In Section 3.4, the curvature of the universe, Qk is allowed 
to float and the fraction of massive neutrinos is once again 
held fixed: 



varied 



fixed 



(r^fc, (Jdm, IjJb, T, Tls, As, 9). 

= 0,woE = -l,r = 0). 



(10) 
(11) 



In Section 3.5, the equation of state of the dark energy is 
varied: 



Table 1. The parameter space probed in our analysis. We as- 
sume a flat prior in each case. We do not vary the values of all 
parameters at the same time; the parameter spaces that we con- 
sider are set out in Section 12.21 



Parameter 


Allowed range 


nk 


-0.3 


- 0.3 


<^dm 


0.01 - 


0.99 


<^ 


0.005 


- 0.1 


Su 


- 


0.5 




-2. 


- 


r 


- 


0.8 


ris 


0.5 - 


1.5 


logio(lOiOAs) 


2.7- 


4.0 


r 


- 


1 


e 


0.5 - 


- 10 



Pviild"" = (t^dm,C^b,WDE,T,ns, As,9). 

In this case, the fixed parameters are: 
Pfixcd°" = (f^K = 0,/. = 0,r = 0). 



(12) 



(13) 



Finally, in Section 3.6, the constraints on tensor modes are 
investigated: 



Pvaricd = (t^dm, t^b, T, Tls, As, T, 9), 

with the fixed parameters given by: 

Pflxed = (f^k=0,/.=0,WDE = -1). 



(14) 



(15) 



Table 1 summarizes the ranges considered for different cos- 
mological parameters when their values are allowed to vary. 



2.3 Constraining parameters 

The prohibitive computational cost of generating CMB 
power spectra and matter transfer functions for all the grid 
points in a multidimensional parameter space has driven 
the development of codes that sample the space selectively, 
guided by the shape of the likelihood surface. We use a 
Markov Chain Monte Carlo (MCMC) approach to search 
the parameter space of the cosmological model (for a recent 
example of the application of the MCMC algorithm to cos- 
mological applications, see Percival 2004) . In brief, this algo- 
rithm involves conducting a series of searches of parameter 
space called chains. The chains are started at widely sepa- 
rated locations within the space. The next link in a chain 
is made in a randomly chosen direction in the parameter 
space. The new link becomes part of the chain if it passes a 
test devised by Metropolis et al. (1953); in summary, links 
for which the likelihood increases are always retained, oth- 
erwise acceptance occurs with a probability that is the ratio 
of likelihoods between the new and old links. If a link is 
rejected, a new randomly generated step is taken in the pa- 
rameter space. This rate of hopping between pairs of points 
in parameter space satisfies the principle of detailed balance, 
so that the chains should asymptotically take up a station- 
ary probability distribution that follows the likelihood sur- 
face. The advantage of this method is that marginalization 
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(i.e. integration of the posterior distribution over uninter- 
esting parameters) is extremely easy: one simply adds up 
the number of links that fall within binned intervals of the 
interesting parameter values (see the appendices in Lewis & 
Bridle 2002). 

The results presented in this paper were generated with 
the publicly available CosmoMC code of Lewis & Bridle 
(2002). We have compared the parameter constraints ob- 
tained with this code with those from an independent code 
written by one of us (WJP), and find excellent agreement 
between the two sets of results. CosmoMC uses the CAMB 
package to compute power spectra for the CMB and matter 
fluctuations (Lewis, Challinor & Lasenby 2000). Our analy- 
sis was carried out in parallel on the Cosmology Machine at 
Durham University. For each parameter set considered, we 
ran twenty separate chains using the Message Passing In- 
terface (MPI) convergence criterion to stop the chains when 
the Gelman and Rubin (1992) statistic R < 1.02, which is a 
significantly more stringent criterion than is usually adopted 
(Verde et al. 2003; Seljak et al. 2005). The length of chain 
generated before the above convergence criterion is achieved 
depends upon the datasets used. For CMB data alone, the 
chains typically have on the order of 10,000 links; in the case 
of CMB plus the 2dFGRS P{k), convergence can be reached 
more quickly. In total, our calculations have accounted for 
the equivalent of more than 30 CPU years on a single pro- 
cessor. 



3 RESULTS 

In this Section, we carry out a systematic study of the con- 
straints placed on the values of cosmological parameters by 
the CMB and LSS datasets listed in Section 12.11 We vary 
three aspects of the comparison: (i) The datasets used. We 
compare constraints obtained from the CMB data alone (Ta- 
ble 2) with those obtained from the CMB data in combi- 
nation with the 2dFGRS power spectrum (Table 3). This 
allows us to see which parameters are constrained more 
strongly when the CMB data are combined with a mea- 
surement of the galaxy power spectrum, (ii) The number of 
parameters varied. We consider models in which 5, 6 or 7 
parameters are allowed to float whilst the other parameters 
are held at fixed values (see Section f2.2l for the definition of 
our parameter spaces), (iii) The combination of parameters. 
In our 7 parameter models, we add one additional param- 
eter to our basic-six set (Eq. ^ and explore how different 
choices for this additional parameter can affect the parame- 
ter constraints. 

Our results are summarized in Tables 2 and 3. In the 
top half of each table, we show the values of the fundamen- 
tal parameters. These are either the range of values derived 
by comparison with a particular dataset or the value that 
a parameter is fixed at in the analysis, as explained in Sec- 
tion In the lower part of the tables, we quote the values 
of other useful parameters (as listed in Eq.|^. These param- 
eters are not varied directly in our analysis. However, their 
values can be derived from the results in the upper half of 
the table, as explained in Section 2.2. 

In Section 3.1, we present the results for a minimal cos- 
mological model with five parameters, the basic-five set. In 
Section 3.2, we consider six parameters, the basic-six set. 



allowing the spectral index of scalar fluctuations to float. 
Sections 3.3 to 3.6 are devoted to seven parameter models, 
with different choices for the 'final' parameter that augments 
the basic-six set, as follows: Section 3.3 the mass fraction 
of massive neutrinos, f^, Section 3.4 non-flat models. Sec- 
tion 3.5 the dark energy equation of state, uide, Section 3.6 
the addition of tensor perturbations. 

In the results tables, unless otherwise stated, we quote 
errors that enclose 68% of the probability around the mean 
value of each parameter. In the subsequent figures showing 
the marginalized posterior likelihood surface for two param- 
eters, the contours mark the locus where — 2A In (£/£max) = 
2.30 and 6.17, corresponding to the 68% and 95% limits re- 
spectively; for the case of a Gaussian likelihood, these con- 
tours correspond to the 'la' and '2a' limits for two degrees 
of freedom. 



3.1 The simplest case — five parameters 

We first concentrate on the simplest possible model that 
gives an accurate description of the data sets, the basic- 
five parameter space defined by Eqs. |1| and |S| This model 
does a remarkably good job of reproducing the CMB data, 
with tight constraints obtained on the values of the sub- 
set of five cosmological parameters varied, as shown by the 
dashed lines in Fig. Qand column 2 of Table 2. It is clear 
from Fig. [Hand Table 3, that, when the 2dFGRS P(k) is 
included, the results show an impressive consistency with 
those obtained from the CMB data alone. For example, in 
the case of the physical density of dark matter, Wdm, the cen- 
tral values derived when comparing to CMB data alone and 
to CMB plus 2dFGRS agree well within the uncertainties. 
However, in a number of cases, there is a significant im- 
provement in the parameter constraints obtained when the 
2dFGRS P{k) data are included. For example, the range of 
Wdm values derived is narrower by a factor of 2 when the 
2dFGRS P{k) is included in the fit, as the LSS data breaks 
the horizon-angle degeneracy arising from CMB models with 
the same position of the first peak in the angular power spec- 
trum (e.g. Percival et al. 2002) . A similar reduction in uncer- 
tainty occurs for the derived parameters as and h. The CMB 
power spectrum is sensitive to the parameter tJdm = ftdmh^ 
whereas the matter P{k) depends on the parameter combi- 
nation Qdmh. The incorporation of P{k) into the analysis 
helps to break the degeneracy between Qdm and h present 
in the theoretical predictions for the CMB, thus tightening 
the constraints on these parameters, as well as on codm- Cole 
et al. (2005) used the 2dFGRS P{k) to place constraints on 
the parameter combinations Q^nh and fib/f^m, and, in con- 
junction with the WMAP temperature power spectrum, on 
Qm- The model that Cole et al. considered is a restricted 
version of our basic-five model (they assumed h = 0.72). It 
is reassuring to note that our results are in excellent agree- 
ment with those obtained by Cole et al.; in particular, we 
confirm their finding of a matter density significantly below 
the canonical ACDM value of flm = 0.3. The success of this 
simple model in describing the current CMB and LSS data 
is remarkable. This 'minimalist model' does a perfectly good 
job of accounting for the form of the most precise probes of 
the cosmological world model that are available to us today. 
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Table 2. Marginalized 68% interval constraints (unless stated otherwise) on cosmological parameters obtained using CMB information 
only for the different hypothesis and parameter sets analysed. The models are defined in Section 2.2. 



b5 b6 b6 + b6 + Qfc b6 + «;de b6 + r 
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Figure 1. Marginalized posterior likelihoods for the cosmological parameters in the basic-five model determined from CMB information 
only (dashed lines) and CMB-|-2dFGRS P{k) (solid lines). The diagonal shows the likelihood for individual parameters; the other panels 
show the likelihood contours for pairs of parameters, marginalizing over the other parameters. The contours show — 2A ln(L/Lniax) = 2.3 
and 6.17. 



3.2 Six parameters — including the scalar spectral 
index 

We now expand our model to allow variations in the scalar 
spectral index, ris, which we call the 'basic-six' parameter 
space (defined by Eqs.|S|and[7|l. Fig.|5|shows the marginal- 
ized likelihoods for this parameter set (along the diagonal), 
together with the two dimensional likelihood contours for 
different combinations of parameters. The results are shown 
using the CMB data alone (dashed lines) and for CMB plus 
2dFGRS P{k) (solid lines). The additional degree of freedom 
gives rise to a well known degeneracy that involves all six pa- 
rameters and which is seen most clearly in the optical depth 
to last scattering, r, and the spectral index and amplitude 
of scalar fiuctuations, ris and As respectively. This degener- 



acy leads to the production of similar power spectra as the 
parameter values, with the exception of Wdm, are increased 
(see Tegmark et al. 2004b for a full description of how the 
degeneracy works in practice). Table 2 shows that in the 
case of the CMB data alone, the results for the best fitting 
parameters in the basic-six case are, for the most part, very 
similar to those obtained for the basic-five parameter set. 
The two exceptions are r and As, for which slightly lower 
values are obtained in the basic-six case. This is also a con- 
sequence of the above degeneracy, since, as the data prefer 
Us < 1, the best fitting values for r and As also decrease. 
Another consequence of the degeneracy is to broaden the 
allowed regions compared with those obtained for the basic- 
five parameter set. The 2dFGRS power spectrum helps to 
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Figure 2. Marginalized posterior likelihoods for the cosmological parameters in the basic-six model determined from CMB information 
only (dashed lines) and CMB plus 2dFGRS P{k) (solid lines). 



break this degeneracy, particularly by tightening the con- 
straints on uidm- The results listed in column 3 of Table 3 
show that the marginalized constraints obtained in this case 
are in complete agreement with those in the CMB only case, 
but with tighter allowed ranges. This reinforces the consis- 
tency of the results obtained from CMB alone and CMB 
plus 2dFGRS P{k) that we found in the basic-five case. 

One particularly remarkable result is the recovered 
value of the spectral index of scalar perturbations, n-s. In the 
case of CMB data alone, we obtain = 0.970^q pggjp ggj) 
where the errors correspond to 68% (95%), fully consis- 
tent with ris = 1. However, with the smaller errors afforded 
by combining the CMB data with the 2dFGRS P{k) data, 
we obtain ria = 0-954^g °23{o o4oj ■ This measurement of the 
scalar perturbation spectral index is consistent with scale 
invariant value ris = 1 at the 95% level. Any detection of 



a deviation from scale invariance would have strong impli- 
cations for the inflationary paradigm, and we discuss this 
result in more detail below in Section 5.2. 

3.3 Six parameters plus the mass fraction of 
massive neutrinos 

Massive neutrinos were ruled out a generation ago as the 
sole constituent of the dark matter, on the basis of N-body 
simulations of the formation of large scale structure in hot 
dark matter universes (Frenk, White & Davis 1983). How- 
ever, interest in massive neutrinos has been resurrected re- 
cently with the resolution of the solar neutrino problem and 
the advent of precision measurements of the galaxy power 
spectrum. The detection of other flavours of neutrino in 
addition to the electron neutrino in the flux of neutrinos 
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Figure 3. The marginalized posterior likelihood in the — 
'^dm plane for the basic-six+/„ parameter set. The dashed lines 
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Figure 4. The marginalized posterior likelihood in the — 
Cue plane for the basic-six+/y parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours show the parameter constraints obtained for 
combined CMB and 2dFGRS P(k) datasets. 



from the Sun suggests that neutrinos can oscillate between 
flavours (Ahmad et al. 2001). This in turn implies that the 
three known types of neutrino have a non-zero mass, al- 
though measurements of the degree of flavour mixing set 
limits on the mass-squared differences between the neutrino 
flavours rather than on their absolute masses. The most ex- 
treme (and perhaps most plausible) case is where the lightest 
mass eigenvalue is negligibly small, in which case the sum 
of neutrino masses is dominated by the heaviest eigenvalue: 
rriv — ~ 0.045 eV (for a recent review see Barger et al. 
2003). The only way in which ^ rrii, can greatly exceed this 
figure is if the mass hierarchy is almost degenerate; we there- 
fore assume three species of equal mass in what follows. Ab- 
solute measurements of neutrino mass can be obtained from 
Tritium beta decay experiments. At present, such experi- 
ments provide a limit on the sum of the neutrino masses of 
J^m,. < 6.6 eV at the 2a level (Weinheimer 2002). 

Currently, the most competitive limits on neutrino 
masses are obtained through the comparison of CMB and 
LSS data with theoretical models (Hu et al. 1998; Elgaroy 
et al. 2002; Hannestad 2002). In the early universe, when 
neutrinos were still relativistic, they free-streamed out of 
density perturbations, damping overdensities in the baryons 
and cold dark matter. This smearing effect stops once neu- 
trinos become non-relativistic; in this case, free-streaming 
only suppresses power on scales smaller than the horizon at 
this epoch, which depends on neutrino mass. 

The CMB temperature power spectrum is only weakly 
dependent on the neutrino mass fraction, f^, since at the 
epoch of last scattering neutrinos with eV masses behave in 
a similar fashion to cold dark matter. So, CMB data alone do 
a poor job of constraining the neutrino mass fraction. More- 
over, the response of the CMB power spectrum to variations 



in is limited to the higher multipoles {£ > 700) and so the 
first year WMAP data alone cannot give good constraints on 
this quantity (see the results of Tegmark et al. 2004b) . Our 
constraints in the CMB only case arise mainly due to data 
other than WMAP which probe smaller angular scales and 
therefore higher multipoles. On the other hand, the impact 
of massive neutrinos on the shape of the matter power spec- 
trum is much more pronounced. The combination of CMB 
data with a measurement of the mass power spectrum can 
therefore give a much tighter constraint on the mass frac- 
tion of neutrinos; the shape of P{k) constrains the value of 
while the CMB data sets the values of the parameters 
that are degenerate with fi,. 

Using CMB data only, we find < 0.182 at 95%. 
When the 2dFGRS P{k) is included this becomes < 
0.105 at 95%. Our results can be converted into constraints 
on the sum of the three neutrino masses using ^ — 
<jJdm/i^ 94.4 eV (assuming standard freezeout and that neu- 
trinos are Majorana particles) to obtain the following lim- 
its: < 2.09 eV at 95% in the CMB only case, 
< 1.16 eV at 95% for CMB data plus the 2dFGRS 

P(fc). 

Elgaroy et al. (2002) used the Percival et al. (2001) mea- 
surement of the 2dFGRS power spectrum to constrain the 
neutrino mass and found '^m.i, < 2.2 eV (95%), assum- 
ing Us = 1 and a restrictive prior on ilm- Our results also 
represent a substantial improvement over those reported by 
Tegmark et al. (2004b), who combined the first year WMAP 
data with the SDSS power spectrum to constrain a similar 
set of parameters to those we consider and found a 95% limit 
of ^ rrii, < 1.7 eV. Our results for provide an important 
illustration of the need to augment the WMAP data, which 
is the most accurate available for £ < 600, with measure- 
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ments conducted at higher angular resolution, allowing sig- 
nificant improvements in the constraints attainable on cer- 
tain parameters. 

It is possible to obtain a stronger limit from CMB-I-LSS 
studies if amplitude information is also used: a neutrino frac- 
tion reduces the overall growth rate as well as changing the 
shape of the matter power spectrum. This constraint was 
used in the year-1 WMAP analysis, and was important in 
reaching the tight constraint of ^ < 0.7 eV (Spergel 
et al. 2003; Verde et al. 2003). This analysis required the 
use of the 2dFGRS bispectrum in addition to P{k) (Verde 
et al. 2002; for a determination with the final 2dFGRS see 
Gaztaiiaga et al. 2005); we have preferred not to use this in- 
formation at the present time, since it has not been subject 
to the same degree of detailed simulation as P{k). The limit 
on the neutrino mass can also be tightened if a measurement 
of the linear theory matter power spectrum is available at 
higher wavenumbers than can be probed with the galaxy 
power spectrum. Seljak et al. (2005) used the power spec- 
trum of the Ly-a forest and the SDSS P{k), with a prior on 
the optical depth to last scattering of r < 0.3 (see later), to 
obtain ^ < 0.42 eV. The extraction of the linear theory 
power spectrum of matter fluctuations from the Lyman-a 
forest remains controversial, so we do not address the use 
of this dataset here (Croft et al. 2002; Gnedin & Hamilton 
2002; McDonald et al. 2005). 

The only work to have reported a measurement of a 
non-zero neutrino mass rather than an upper limit is Allen 
et al. (2003). These authors combined galaxy cluster data 
with CMB data and an earlier version of the 2dFGRS P{k) 
measured by Percival et al. (2001). The cluster data used 
by these authors was the gas fraction and the X-ray lumi- 
nosity function; both quantities are much more difficult to 
model than the CMB and LSS data that we consider here. 
Although their results show a stronger signal upon the inclu- 
sion of the galaxy cluster data, there is still the suggestion 
of a non-zero neutrino mass fraction even with the CMB 
and 2dFGRS P{k) data alone, showing that this conclusion 
is not due exclusively to the use of the X-ray data. The pa- 
rameter space explored by Allen et al. differs from the one 
considered in this section, since it includes tensor modes. 
The tensor modes contribute to the low multipole part of 
the CMB spectrum and their inclusion can drive down the 
amplitude of the scalar perturbations on these scales. This 
in turn can lead to an increase in the recovered value of the 
scalar spectral index, ris, with the consequence that fi, in- 
creases to compensate, thus maintaining the power in the 
mass distribution at high k. This degeneracy in the fv — r 
plane produces a higher one-dimensional marginalized con- 
straint on the neutrino mass fraction. 

Fig. |3] and Fig.|3]show the impact of including the 2dF- 
GRS P{k) data on the fv—uidm and /,y — Ode constraints. In 
the CMB only case, the incorporation of fv into the param- 
eter space causes the uncertainty in all the parameters to 
grow. This is particularly noticeable for SIdb, for which the 
errors are twice as big as they were for the basic-six param- 
eter set with fv = 0. When the 2dFGRS power spectrum is 
added to the analysis, the allowed ranges of these parameters 
are dramatically reduced, with particularly tight constraints 
resulting on ojdm and JIde; this clearly demonstrates the im- 
portance of including LSS data to obtain precise constraints 
on these parameters. 
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Figure 5. The one dimensional marginalized posterior likeliliood 
for Qfe for CMB data only (dashed line), CMB plus 2dFGRS P{k) 
(solid line), and CMB plus 2dFGRS P{k), with a prior on the 
optical depth of t < 0.3 (dot-dashed line). Closed models have 
Cfe <0. 



3.4 Six parameters plus the curvature of the 
universe: non-flat models 

There is a strong theoretical prejudice that we live in a fiat 
universe with Q.k ~ 0. The first detections of the acoustic 
peaks in the CMB temperature power spectrum, the loca- 
tion of which is a measure of the geometry of the universe, 
showed that the universe is close to being fiat (de Bernardis 
et al. 2000). These results served to reinforce the prejudice 
that the curvature of the universe must be exactly zero - 
and it is true that, to date, no work has found any strong 
indication of a significant deviation from Qk = 0. However, 
as the flatness of the universe is one of the most important 
predictions of inflationary models, this assumption must be 
properly tested against new datasets. We must bear in mind, 
when comparing values reported for cosmological parame- 
ters, that many works simply assume Qk ~ 0. Other param- 
eters, for example the scalar spectral index, are sensitive to 
the prior assumed for Qk- 

We plot the marginalized likelihood function for ilk iu 
Fig. 13 for different datasets. The dashed curve shows the 
results for the CMB data alone, reminding us of the well- 
known (but frequently forgotten) result that the CMB data 
alone do not require a flat universe. Even though values of 
Qk > (open models) are practically ruled out, a wide range 
of closed models is allowed, with the best fitting value given 
by Qk = -OmA+llll at 68% (95%) confidence. The 

solid line in Fig. |3 shows how incorporating the 2dFGRS 
power spectrum helps to tighten the constraints on Qk - The 
addition of power spectrum information helps to break the 
geometrical degeneracy between nm and Qde (see Fig. |3 
and the final paragraph of this subsection). This is one of 
the most important effects of the incorporation of LSS in- 
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Figure 6. The marginalized posterior likelihood in the fi^ — r 
plane for the basic-six plus Hj. parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours correspond to the constraints obtained in the 
CMB plus 2dFGRS P(fe) case. 



Figure 7. The marginalized posterior likelihood in the f!m— f^DE 
plane for the basic-six plus fl^ parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours correspond to constraints in the CMB plus 
2dFGRS P{k) case. 



formation into the analysis. In the CMB plus 2dFGRS P{k) 

i n n r,oo-l-0.018 (0.032) 

case, we get fit = -0.029_(, (,jg pjgj . 

It is particularly important to note the effect that the 
prior on the optical depth to the last scattering surface, r, 
has on the inferred value of the curvature of the universe. 
Fig. |S| shows the constraints in the Qk — t plane. The ad- 
dition of the 2dFGRS power spectrum shrinks the allowed 
region by tightening up the constraints on flk, but the re- 
sulting likelihood contours show a clear degeneracy between 
the two parameters, with the high values of r preferring 
more negative values of fifc. This degeneracy is responsi- 
ble for the broad error bars on these parameters. If one 
adopts a restrictive prior on the optical depth of r < 0.3, 
as recommended by the WMAP team based on the lack of 
a large signal in polarization autocorrelation, then the re- 
sults for Qk are more in line with those in the literature, 
as shown by the dot-dashed line in Fig. 13 In this case we 
find flk ~ ~0.015^g pi[j (0 020) ^'-'^ combined CMB plus 
2dFGRS datasets. We shall return to the issue of the choice 
of prior for the optical depth in Section f4. 41 

Finally, we highlight the constraints on the densities of 
dark matter and dark energy obtained when the assumption 
of a flat universe is dropped. Fig.|2|shows the results for the 
case of CMB data alone (dashed lines) and for CMB data 
plus the 2dFGRS P{k) (solid lines). As we have seen in 
several previous examples, there is a dramatic improvement 
in the quality of the constraints on these parameters once 
the galaxy clustering data is incorporated into the analysis. 
There is compelling evidence for a dark energy component 
in the universe. 



3.5 Six parameters plus the dark energy equation 
of state 

Over the past decade, mounting evidence has been presented 
for the accelerating expansion of the Universe, based on 
the interpretation of the Hubble diagram of Type- la super- 
novae (Perlmutter et al. 1999; Riess et al 2004). Independent 
support for the presence of a dynamically dominant, nega- 
tive pressure component in the energy-density budget of the 
Universe has also come from fitting cosmological models to 
CMB and LSS datasets (see Section r-{.4l and Efstathiou et al. 
2002). Although we can infer the presence of this component, 
dubbed dark energy, we know practically nothing about its 
nature. A plethora of theoretical models have been proposed 
for the dark energy (e.g. see the review by Sahni 2004). One 
of the key properties of the dark energy which can be used 
to pare down the market of possible models is the equation 
of state of the dark energy, that is the ratio of its pressure 
to density, iude- 

Until now we have assumed that the dark energy com- 
ponent corresponds to the cosmological constant, with a 
fixed equation of state specified by woe = —1- However, 
this is only one manifestation of the many possible forms 
that the dark energy could take. Any component with an 
equation of state wde < — 1/3 will result in an accelerating 
rate of expansion today. In this section, we explore dark en- 
ergy models with a constant equation of state, allowing for 
variations in the redshift independent value of woe- We also 
consider models with wde < —1, sometimes referred to as 
'phantom energy'. 

Fig.|H|shows the marginalized constraints in the uidb — 
Qui plane. In the CMB only case, we find uide = — O.QSIq 47, 
consistent with a cosmological constant. When the 2dFGRS 
power spectrum is included in the analysis, the preferred 
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Figure 8. The marginalized posterior likelihood in the wde 
plane for the basic six plus iode parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours show the corresponding constraints obtained 
in the CMB+2dFGRS P(fc) case. 



Figure 9. The marginalized posterior likelihood in the ris — r 
plane for the basic-six plus r parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours show the corresponding results in the CMB 
plus 2dFGRS P{k) case. 



value increases somewhat to tode = — 0.85_q '^y. If we also in- 
clude the supernova type la data from Riess et al. (2004), our 
result scarcely changes, with iude = — 0.87^° ijj- Phantom 
energy models are permitted in the case of CMB data only, 
with the 95% limit on the equation of state of wde > —1.66. 
However, once the 2dFGRS P{k) and supernovae type-la 
data are included, the allowed region shrinks to a smaller 
zone with wde > —1.19 at 95%, showing that phantom en- 
ergy models are disfavoured by the currently available data. 
These results show that the data prefers lower values of wde 
than suggested by previous work using the SDSS power spec- 
trum (MacTavish et al. 2005). Our results are consistent 
with the dark energy taking the form of a cosmological con- 
stant. We will discuss this point further in Section |S| 

3.6 Six parameters plus non-zero tensor modes 

We now add the ratio of the amplitude tensor to scalar per- 
turbations, r, to the basic-six parameter set. This case is an 
important one to consider as tensor modes are predicted 
to be present in many inflationary models. Moreover, as 
we shall see, several cosmological parameters are degenerate 
with r and in the literature tensor modes have often been 
ignored when presenting constraints on these parameters. 

The constraints imposed on r by CMB information 
alone are r < 0.52 at 95%. Including the 2dFGRS P{k) data 
reduces the importance of tensors slightly, yielding r < 0.41 
at 95%. Fig.Elshows the two-dimensional marginalized like- 
lihood contours in the rig — r plane for the cases of CMB 
data only (dashed lines) and CMB plus the 2dFGRS P{k) 
(solid lines). Tensor modes contribute to the CMB temper- 
ature power spectrum only on large angular scales, leading 
to a reduction in the scalar perturbations on these scales to 



match the observations. In order to maintain the amplitude 
of scalar perturbations on smaller angular scales, an increase 
in the scalar spectral index, ris, is required. This degeneracy 
results in a broader allowed range for Us than in the case 
where only scalar modes are considered. 

The constraints on r and rig can be translated into the 
horizon flow parameters, ei and £2 , using the relations given 
by Mukhanov et al. (1992): 

1-ne = 2ei+e2 (16) 
r = 16ei. (17) 

The horizon flow parameters are related to derivatives of 
the Hubble parameter during inflation (Schwarz et al. 2001). 
Leach & Liddle (2003a) give equations relating the horizon 
flow parameters to the derivatives of the inflaton potential 
and discuss the motivation for the truncation of the slow- 
roll expansion after £2. The constraints on the horizon flow 
parameters are shown in Fig. IIUI The degeneracy between r 
and Us translates into a degeneracy in ei and e2. 

If we restrict our attention to monomial inflation, i.e. 
potentials of the form V oc (jj" , then the horizon flow pa- 
rameters can be related to the power law index, a, and the 
number of e-folds of inflation for the scale considered. A'', by 
the simple relations (Leach & Liddle 2003b): 

« = ti- (18) 

N . -l). («) 

To obtain constraints on these new parameters, we have 
translated our results for ei and £2 into the a — N plane. 
In doing so, we have restricted our attention to the region 
where e2 > 0, following Liddle & Leach (2003b), who argue 
that this part of the horizon flow parameter space contains 
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Figure 10. The marginalized posterior likelihood in the ei — £2 
plane for the basic-six plus r parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours correspond to the results obtained in the CMB 
plus 2dFGRS P(k) case. 



Figure 11. The marginalized posterior likelihood in the a — N 
plane for the basic-six plus r parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours show the results in the CMB plus 2dFGRS 
P{k) case. 



the most likely models in which inflation will end naturally 
with a violation of the slow-roll approximation. Our results 
are plotted in Fig.[TTl We find that a < 2.33 at 95% for CMB 
data alone and a < 2.27 (95%) for CMB plus the 2dFGRS 
P(k). To obtain this result, we have followed Seljak et al. 
(2005) and take into account the maximum number of e- 
folds, A^max = 60, of slow-roU inflation experienced at the 
pivot scale k = 0.05 Mpc~^, thus further restricting the sec- 
ond horizon flow parameter, £2 > 0.0167. (Note that Leach 
& Liddle 2003b use a different pivot scale to ours.) Our con- 
straint on a implies that the A</>^ inflation model is ruled 
out. This is the first time that the CMB data alone have 
been of sufficient quality to completely reject this model. 
Seljak et al. (2005) reached a similar conclusion using dif- 
ferent datasets: the WMAP data, the SDSS galaxy power 
spectrum and either the power spectrum of the Ly-a forest 
or the amplitude of the matter power spectrum, as inferred 
from the bias of SDSS galaxies. 

We now turn our attention to power-law inflationary 
models, in which the scale factor of the universe grows with 
time as a ccf with p > 1. In such models, the horizon flow 
parameters are simply given by 

ei = i (20) 
P 

ei = i > 2. (21) 

Substitution of these relations in Eo. 1171 gives a relation be- 
tween Tis and r: 

r = 8(1 -Tie). (22) 

To analyse this kind of model, we ran a new set of chains 
fixing the tensor to scalar ratio using Eq. 1221 In this case 
we get Us = 0.978l° olo for the CMB data alone and 



Us = G.'d7&2Zo^lll in the CMB plus 2dFGRS case. The 
constraints on r are also tighter, with r < 0.31 at 95%. We 
note that the best fitting values for the horizon fiow param- 
eters of ei = 0.0123lg;^{J^^ (corresponding to p = Slljf ) 
and 62 = — 0.004lg Q4Q are in complete agreement with the 
power law inflation picture. 



4 THE ROLE OF PRIORS 

It is often claimed that we have entered an era of precision 
cosmology, in which the values of the cosmological param- 
eters are known with high accuracy. The CMB measure- 
ments alone go a long way towards realizing this ideal, but 
ultimately fall short due to the presence of well known de- 
generacies between the cosmological parameters (Efstathiou 
& Bond 1999). Some of these degeneracies can be broken 
with the incorporation of other information into the analysis 
(such as, for example, LSS, SN la or the power spectrum of 
the Ly-Q forest). However, many degeneracies remain even 
after the addition of these datasets. Another way to break 
degeneracies is by imposing priors on parameters, which can 
have implications for the derived parameter constraints (see, 
for example. Bridle et al. 2003a). In this section we revisit 
the constraints obtained for the different parameter sets and 
priors and assess which of our results are the most robust. 



4.1 The baryon density 

One of the most important achievements of modern cosmol- 
ogy is the agreement between the value of the physical den- 
sity of baryons determined from CMB data and that inferred 



14 Ariel G. Sanchez et al. 



from Big Bang Nucleosynthesis (BBN) arguments and dis- 
tant quasar absorption spectra. In the present analysis, we 
obtain a value for the baryon density of ujh = 0.02291q qois 
from the CMB data alone that is consistent with the lat- 
est constraint from BBN: c^b = 0.022 ± 0.002 (Cuoco et al. 
2004). This agreement is reinforced when the 2dFGRS P{k) 
is added to the analysis, with coh ~ 0.02251q qqJ^q . The vari- 
ation in the value of ujh obtained between the different pa- 
rameter sets and priors that we have analysed is smaller 
than the la error bars, showing the robustness of this re- 
sult. This level of agreement is all the more remarkable when 
one considers the quite different epochs to which the various 
datasets relate: BBN is a theory that describes processes oc- 
curring in the very early universe, just a few minutes after 
the Big Bang, while the CMB maps the universe as it was 
a few hundred thousand years after the Big Bang, and the 
galaxy power spectrum refers to the present day universe, 
over 13 billion years later. The fact that we can tell a co- 
herent story over such a huge baseline in time and physical 
conditions provides an impressive verification of our cosmo- 
logical model. 
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4.2 The dark matter density 

A scan across the third rows of Tables 2 and 3 shows that 
the value of ojdm is largely insensitive to the priors applied to 
the other parameters. The one exception is when the flatness 
prior, Qk ~ 0, is relaxed, in which case we obtain a smaller 
value for oJdm with larger errors. The constraints obtained 
on LJdm in the CMB only and CMB plus 2dFGRS P{k) cases 
are fully consistent. 

The implications of the value of uidm for the matter den- 
sity Qui do, however, depend on the priors implemented. In 
the basic-six plus fi, parameter set, for the CMB plus 2dF- 
GRS case we find S7m = 0.282 ± 0.040, but it can be as low 
as Qm = 0.224 ±0.022 for the basic-six plus r parameter set. 
With the exception of the case of non-zero neutrino mass, 
all our estimates of flm lie significantly below the standard 
choice of 0.3. Fig. 1121 illustrates how the choice of param- 
eter space affects the results obtained. The constraints in 
the ilm — o"8 plane in the basic-six parameter set are tighter 
than those obtained when the neutrino fraction, is incor- 
porated into the analysis; for the latter case, a bigger region 
with lower values of as is allowed by the data. A similar 
situation can be seen in Fig. 1131 for the Sim — h plane. The 
values of h preferred by the data are lower when non-flat 
models are considered in the analysis. These discrepancies 
cause differences in the marginalized results obtained for 
these parameters. This situation occurs in many other cases 
and in general the influence of the parameter set is non- 
negligible. For this reason, constraints on a given parameter 
should always be quoted together with the parameter space 
explored in the analysis. Fig. 1131 shows the la limits on the 
Hubble constant derived by the HST Key Project (Freed- 
man et al. 2001). These constraints on h are substantially 
broader than those obtained from the CMB plus 2dFGRS 
power spectrum, showing that including the Key Project 
measurement as a prior would have little impact on our re- 
sults. 



Figure 12. The marginalized posterior likelihood in the Qrn — o"8 
plane obtained using CMB plus 2dFGRS information for different 
parameter sets. The solid lines correspond to the 68% and 95% 
contours obtained for the basic-six parameter set. The dashed 
lines correspond to the results obtained when the neutrino frac- 
tion fi, is also allowed to vary. The dot-dashed lines show con- 
straints from weak lensing measurements from Hoekstra et al. 
(2002). 
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Figure 13. The marginalized posterior likelihood in the Qm — h 
plane obtained using CMB plus 2dFGRS information for different 
parameter sets. The solid lines correspond to the 68% and 95% 
contours obtained for the basic-six parameter set. The dashed 
lines show the results obtained when non-fiat models are consid- 
ered (basic-six plus Oj.). The dot-dashed lines show the la con- 
straint on h from the HST Key Project (Freedman et al. 2001). 
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Figure 14. The marginalized posterior likelihood in the fi^ — to 
plane for the basic-six plus Hj. parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours correspond to the constraints obtained in the 
CMB plus 2dFGRS case. 



Figure 15. The marginalized posterior likelihood in the fi^ — ris 
plane for the basic six plus Hj. parameter set. The dashed lines 
show the 68% and 95% contours obtained in the CMB only case. 
The solid contours correspond to the constraints obtained in the 
CMB plus 2dFGRS case. 



4.3 The amplitude of fluctuations 

When constraining the values of the cosmological parame- 
ters, we only use information from the shape of the galaxy 
power spectrum and not from its amplitude. Therefore the 
constraints on the amplitude of density fluctuations come 
principally from the CMB data, with the LSS data playing 
an indirect role by tightening the constraints on parame- 
ters which yield degenerate predictions for the CMB. The 
recovered values of cts range from erg = 0.678lg g!^^ for the 
basic-six plus parameter space to erg = CSlTlg gi^g for 
non-flat models (basic-six plus fifc). Adding more data sets, 
such as the X-ray cluster luminosity function, or other mea- 
surements of the amplitude of fluctuations may help to im- 
prove the constraint on as, but the theoretical modelling 
of these observational datasets is less straightforward. In 
Fig. 1121 we compare our results with constraints from mea- 
surements of weak lensing from Hoekstra et al. (2002) . Henry 
(2004) used the temperature function of X-ray clusters to 
find as = 0.66 ± 0.16, in good agreement with the basic-six 
plus neutrino mass fraction model. Similar constraints have 
been obtained by other groups (Bacon et al. 2003; Heymans 
et al. 2005). 

4.4 The optical depth 

The optical depth to the last scattering surface has an 
important effect on nearly all other cosmological parame- 
ters. The signal for r > comes from the temperature- 
polarization cross-correlation on large angular scales. In- 
triguingly, Hansen et al. (2004) performed a temperature- 
polarization analysis of the first year WMAP data for the 
northern and southern hemispheres separately and found 
that, whereas the northern hemisphere points to r = 0, the 



southern hemisphere prefers a value of r = 0.24_q q7, incon- 
sistent with T = at the 2cr level, with the suggestion that 
the signal for r > may be due to foreground structures in 
the southern hemisphere. 

In their analysis of the first year WMAP data Spergel 
et al. (2003) imposed a prior of r < 0.3, justifying this by 
the need to avoid 'unphysical' regions of parameter space. 
In Section 3.4 we demonstrated, as previously pointed out 
by Tegmark et al. (2004b), the strong effect this prior has 
on our results. In particular, the r < 0.3 prior is required to 
reconcile the constraints on Q,k with the flatness prediction 
from inflation in the basic-six plus Q.^ parameter set, and to 
produce tight constraints on neutrino masses in the basic-six 
plus fv case. The situation should improve with the release 
of the second and subsequent years of data from WMAP, 
which will be able to produce improved polarization maps. 
In the meantime, the effect of this important prior must 
be borne in mind when interpreting the results from multi- 
parameter analyses. 

4.5 The flatness prior 

The prior of fJ^ = is widely used when constraining cos- 
mological parameters. It is important to remember that, if 
the assumption of ffatness is relaxed, the preferred value is 
actually fli, < and only marginally consistent with fifc = 0. 
The assumption of flatness has a major impact in the values 
of many cosmological parameters. 

The value obtained for the age of the universe, to, shows 
an important change when Q.k is allowed to float, and is only 
marginally consistent with the values found for fi^ =0. 
Fig. 1141 shows the marginalized two-dimensional likelihood 
contours in the Q,k — io plane. There is a clear degeneracy 
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between these two parameters, with lower values of Qk pre- 
ferring higher values of to; the incorporation of the 2dFGRS 
P{k) data does not break this degeneracy completely. The 
same degeneracy can be seen in the Q.k — h plane, which im- 
phes that a prior on the Hubble constant from the HST key 
project (Freedman et al. 2001) may improve the situation, 
but even then the constraints on these parameters are less 
robust than in the fiat case. 

The scalar spectral index, ns, also merits special atten- 
tion. In Section 3.2, we pointed out that the constraint on 
ris in the basic-six model is only marginally consistent with 
ris = 1; this spectrum is formally excluded at the 95% confi- 
dence level. Fig. lldl shows the marginalized two-dimensional 
likelihood contours in the fifc — ns plane. When CMB in- 
formation alone is used, there is a wide allowed region that 
shrinks considerably when the 2dFGRS power spectrum is 
included. In the latter case, there is a correlation between 
the parameters which makes the constraints on ris much 
broader than those obtained for the special case of fife = 0. 
Taking into account that the evidence for ris < 1 is weaker 
once more general parameter sets are considered (such as, 
for example, the basic-six plus fv set) and that the current 
data has a slight preference for closed models (even when 
the prior of r < 0.3 is applied), we advocate caution before 
claiming a detection of a significant deviation from scale in- 
variance. 

4.6 Tensor modes 

Another commonly applied prior is the assumption that ten- 
sor modes can be neglected. It is important to include the 
amplitude of tensor modes as a free parameter, not only 
because this has strong implications for infiationary mod- 
els, but also because many other parameters are degenerate 
with the amplitude of tensors, resulting in the growth of the 
allowed regions for these parameters. The parameters that 
are most strongly influenced by the assumption about tensor 
modes are ris, Om and h. 



5 BEYOND THE SIMPLEST MODEL 

5.1 How many parameters should float? 

We have shown that a model in which five parameters are al- 
lowed to vary gives a good description of the CMB and LSS 
datasets. We then went on to explore six and seven parame- 
ter sets, finding that, in some cases, the results obtained for 
certain parameters depended upon the choice of parameters 
varied. But are we justified in adding extra free parameters 
to our basic-five set? 

The simplest way to make an objective assessment of 
different models is to establish whether or not they afford a 
better description of the data, which is usually done by com- 
puting likelihood ratios. However, it is important to compen- 
sate for the fact that adding extra parameters should nec- 
essarily improve the fit to the data. Liddle (2004) has advo- 
cated the use of two simple statistics that quantify the level 
of improvement in the description of cosmological datasets 
as new free parameters are added to the theoretical models. 
These statistics allow us to ascertain whether the addition 
of a new parameter is justified, i.e. does it produce a better 



Table 4. The number of parameters, the likelihood and the 
values of the AIC and BIC statistics for the various models analysed 
in this paper. In all cases, A^jata = 1403. 
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than expected enhancement in the accuracy with which the 
data is reproduced? The statistics, the Akaike Iifformation 
Criterion (AIC; Akaike 1974) and the Bayesian Information 
Criterion (BIC; Schwarz 1978) have a long track record of ap- 
plication in other branches of physics, but have largely been 
ignored in cosmology. The definitions of the two statistics 
are straightforward: 

AIC = -21n(£) + 2iVpa,., (23) 
BIC = -21n(/:) +Arpa.ln(iVdata), (24) 

where jC is the maximum likelihood, A^par is the number of 
parameters varied in the model and A'data is the number of 
datapoints included in the analysis. The model that best de- 
scribes the data with the most economical use of parameters 
is the one that minimizes these quantities. In both expres- 
sions, the first term favours models which provide better fits 
to the data, while the second penalizes large numbers of pa- 
rameters. We note that, as the value of In(Adata) > 2 in 
our analysis, the BIC actually gives a higher penalty to the 
number of free parameters than is the case for the AIC. 

Table 4 provides a summary of the number of parame- 
ters, the likelihood and the values of the AIC and BIC statis- 
tics for the models considered in this paper. The addition of 
extra free parameters does of course lead to an increase in 
the likelihood of the description of the data by the model. 
The message conveyed by the value of the AIC statistic is less 
clear. All models show a slight decrease in the value of the 
AIC statistic compared with the basic-five set, but the BIC 
statistic paints a quite different picture. Liddle reports that 
a difference in the BIC of 2 should be regarded as 'positive 
evidence' and of 6 as 'strong evidence' in favour of the model 
with the smaller value of BIC. Therefore, there is apparently 
'positive evidence' that we should not expand the basic-five 
model to allow the scalar spectral index to fioat and 'strong 
evidence' that we really should not have burdened the reader 
with the basic-six plus one more free parameter models. 

This conclusion, and indeed the basic BIC formula it- 
self, appears to disagree with the approach of this paper. 
For data with Gaussian errors, the addition of a further pa- 
rameter would be expected to reduce —2\n£ by one via 
the usual 'degree of freedom' rule (although note that this 
strictly applies only to parameters linearly related to the 
data, such as polynomial expansion coefficients). Further- 
more, the reduction in — 21n£ should be distributed as 
with 1 degree of freedom if the new parameter is not part of 
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Figure 16. Model fits to the CMB datasets (top panels) and 2dFGRS P{k) (bottom panels) in the cases of the basic-five (solid line) 
and basic-six (dashed line) models, with the best fitting parameter values listed in Table 3. The model P{k) have been convolved with 
the window function of the 2dFGRS from Cole et al. (2005). In the right-hand panels, the model curves and datapoints have been divided 
by the best fitting basic-five model to expand the y-axis. 



the true model, independent of the number of data points. 
A reduction in — 21n£ of 4 therefore amounts to rejection 
at 95% confidence of the hypothesis that a new parameter 
is not required. Thus, the fact that allowing deviations from 
scale invariance reduces — 2 In £ by 3.7 amounts to marginal 
evidence for the reality of tilt. This reasoning matches the 
AIC approach quite well, as long as the coefficient 2 in the 
2A'par term is regarded as being adjustable according to the 
significance level of interest. 

The BIG statistic is an approximate form of the 
'Bayesian evidence' (Hobson, Bridle & Lahav 2002, Liddle 
2004, Trotta 2005). One of the conditions that must be sat- 



isfied in order for the BIG to be a good approximation to the 
Bayesian evidence is the independence of the data points un- 
der consideration. By setting A'data = 1403 in the definition 
of the BIG in Eg. 1241 we are effectively treating all of the data 
points used in our analysis as independent. Our calculation 
of the BIG therefore gives an overly pessimistic impression 
of the impact of adding of new parameters. If, for example, 
an eigenmode analysis or 'radical' data compression tech- 
nique was applied to the full set of CMB data points used 
in our analysis, this would produce a much smaller set of 
genuinely independent data points which fully describe the 
CMB measurements (Bond, Jaffe & Knox 2000). The values 
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of the AIC and BIG statistics would become closer if recom- 
puted for this 'reduced' set of datapoints. One might argue 
that one should compute the Bayesian evidence rather than 
approximations such as the BIG. There are two reasons why 
we have not done this. Firstly, the Bayesian evidence is hard 
to compute accurately using MCMC techniques, although 
fast algorithms are under development (Mukherjee, Parkin- 
son & Liddle 2005). Secondly, the definition of the prior on 
a parameter is part of the model tested in the Bayesian ev- 
idence approach and we believe that this is a weak point 
in the method for the following reason. Since the choice of 
prior is arbitrary to some extent, it is possible in principle 
to select a prior such that the Bayesian evidence increases 
upon the addition of new parameters. We prefer the effec- 
tively frequentist argument of simply requiring a reduction 
in — 21n£ of order unity in order to claim the detection of 
another degree of freedom. 
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5.2 Details of the evidence for tilt 

Putting aside the caveat raised by the BIG statistic, it is im- 
portant to look at the basic-five and basic-six model results 
for the scalar spectral index in more detail, as they have 
important implications for infiationary models. To recap, in 
Section 3.2 we set Us = 1, i.e. the scale invariant value of 
the spectral index for primordial scalar fiuctuations. In Sec- 
tion 3.3, we treated the spectral index as a free parameter 
and found that Us = 1 was on the 95% limit. Fig. 1161 shows 
the best fitting models to the CMB temperature power spec- 
trum data (upper panels) and the 2dFGRS P{k) (lower pan- 
els) for the basic-five (solid lines) and basic-six (dashed lines) 
models. The difference between the two models is small and 
comes mostly from scales beyond those probed by WMAP; 
this is quantified in the right-hand panels in which the mod- 
els and datapoints have been divided by the ris = 1 model. 
Fig. 1171 shows the likelihood quotients between the basic- 
five and basic-six models for each dataset separately. It is 
clear that the datasets primarily responsible for driving the 
scalar spectral index away from the scale invariant value are 
the CBI measurements and the 2dFGRS P{k): the basic- 
six model represents only a modest improvement over the 
basic-five model in its description of the WMAP and VSA 
datasets, while the addition of an extra parameter makes 
very little difference to how well the ACBAR results are re- 
produced. Nevertheless, there is an impressive consistency 
between the various datasets: a systematic error in a sin- 
gle one of these might have resulted in an improved overall 
likelihood on the introduction of tilt, but at the price of a 
poorer fit to some of the correct datasets. This is not what 
we see: addition of the 2dFGRS strengthens a weak trend 
already present in the CMB data. Even so, the overall re- 
sult remains tantalisingly placed in terms of its statistical 
significance: 95% confidence is not sufficient to claim firm 
detection of an effect of this importance. The best that can 
be said is that even a modest amount of extra data could 
easily move things into the territory of firm detection. The 
largest predicted deviations from Us = 1 occur around the 
third CMB peak, at ^ ~ 800, and the data here may be 
expected to improve rapidly. 



Figure 17. The likelihood ratios of the basic-five model to the 
basic-six model plotted in Fig |16l for the individual datasets used 
in our analysis. 

6 COMPARISON WITH CONSTRAINTS 

OBTAINED USING THE CMB DATA AND 
THE SDSS POWER SPECTRUM 

In this section, we replace the 2dFGRS P{k) measured by 
Cole et al. (2005) with the power spectrum of SDSS galaxies 
estimated by Tegmark et al. (2004a) and examine the impact 
that this change has upon the values of the recovered cos- 
mological parameters. There are a number of differences be- 
tween the two measurements of the galaxy power spectrum. 
Firstly, the SDSS is a red-selected survey, while the 2dFGRS 
is blue-selected. Secondly, Tegmark et al. used a sophisti- 
cated eigenmode deconvolution apparatus to attempt to re- 
move the effects of the survey window and redshift-space 
distortions; in contrast. Cole et al. used a simpler Fourier 
approach that compares to window-convolved models and 
quantifies redshift-space effects directly by comparison with 
realistic simulations. 

Tegmark et al. (2004b) used the WMAP first year data 
and the SDSS galaxy power spectrum to constrain cos- 
mological parameters. These authors modelled the galaxy 
power spectrum with a non-linear model for the matter 
fiuctuations multiplied by a scale independent bias factor. 
The power spectrum data were used on scales larger than 
k < 0.20/iMpc~^. It is not clear that a constant bias is 
a good approximation on scales for which the density fluc- 
tuations have become non-linear. We adopt a simpler ap- 
proach and assume that the galaxy power spectrum can 
be related to the linear perturbation theory power spec- 
trum of the mass through a constant shift in amplitude. 
As discussed earlier, the simulations used by Cole et al. 
indicate that redshift-space effects and other nonlinearities 
are unimportant for the 2dFGRS to our imposed limit of 
fcmax = 0.15/iMpc"^ 

We repeat the study of parameter space previously car- 
ried out using the CMB data plus the 2dFGRS P{k) and 




Figure 18. The marginalized one-dimensional posterior likelihood in the basic-six parameter space obtained for CMB information only 
(dashed Unes), CMB plus the 2dFGRS P{k) (soUd lines) and CMB plus the SDSS (dot-dashed lines). 



present our results using the SDSS P{k) instead in Table 5. 
For the most part, the results obtained with the SDSS P{k) 
are compatible with those found using the 2dFGRS P{k). 
There are, however, some cases in which the results are quite 
different. This point is illustrated using the results of the 
basic-six model in Fig. 1181 In this plot, we compare the pa- 
rameter constraints obtained using CMB data alone (dashed 
line) with the resuHs for CMB data plus the 2dFGRS P{k) 
(solid line) and for CMB data plus the SDSS P{k) (dot- 
dashed line). For three out of the six parameters presented, 
ujh , T and Us , there is impressive agreement between the sets 
of results; the peaks in the likelihood distributions coincide 
in the CMB only and CMB plus P{k) cases, and the results 
are consistent to better than the 68% confidence intervals. 
The agreement between the sets of results for the scalar spec- 
tral index in particular is excellent. However, for the other 
three parameters plotted, Wdm, h and as, the agreement is 
less impressive; the differences in the recovered values of h 
and as are driven by the change in uidm- The peak in the 
likelihood distributions for the CMB only and CMB plus 



2dFGRS P{k) cases are in good agreement, as remarked 
upon in Section 3. There is a clear discrepancy, however, 
with the preferred parameter values when using CMB data 
plus the SDSS P{k). This is most marked for the physical 
density of dark matter, uJd-m- Cole et al. (2005) noted that 
the SDSS P(fc) has a slightly bluer slope than that of the 
2dFGRS, favouring higher values of 57ni. 

There are two other notable discrepancies between the 
resuhs obtained with the SDSS and 2dFGRS P{k) in our 
basic-six plus one additional free parameter models. When 
the assumption of a flat universe is relaxed, we find that 
the constraints on $1^ are weaker in the SDSS case, Q.^ = 
— 0.070 [l^ggggjQQ^gj ; the allowcd range is nearly twice as 
broad as in the case of the 2dFGRS P(k). This is because 
on the scales used in our analysis, the SDSS power spec- 
trum does a poorer job of constraining compared with 
the 2dFGRS P{k), and hence is not as effective at breaking 
the geometrical degeneracy. In the basic-six plus dark energy 
equation of state parameter set, we find tode ~ ^0.451q23 
using the SDSS P{k), much higher than we found in the case 
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Table 5. The marginalized 68% interval constraints (unless otherwise stated) on cosmological parameters obtained using CMB data and 
the SDSS galaxy power spectrum for different parameter sets. 
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of the 2dFGRS P{k) and inconsistent with a cosmological 
constant. If wc also include the SNIa data from Riess et al. 
(2004), then we obtain a value for the equation of state that 
is consistent with our previous results: wde = — O.SQlg'ig. 
Again, the discrepancy in the result for the equation of state 
can be traced back to the preferred values of ujdm- Wc note 
that MacTavish et al. (2005) find similar results to ours for 
the equation of state of the dark energy using the SDSS 
galaxy power spectrum. Fig. 8 shows the degeneracy in the 
f2m — Wde plane for CMB data alone. Adding information 
from the galaxy power spectrum breaks this degeneracy. If 
the galaxy P(fc) data prefer a high value of f2m, as is the 
case for the SDSS data, then a high value of wde will result. 



It may be that these differences between 2dFGRS and 
SDSS amount to no more than an unlucky amount of cos- 
mic variance, but clearly it would be more reassuring if the 
results showed greater consistency. It will therefore be im- 
portant to see each dataset subjected to analysis by a variety 
of algorithms and codes, as happened following the Percival 
et al. (2001) 2dFGRS power-spectrum analysis (Tegniark, 
Zaldarriaga & Hamilton 2001). This older comparison found 
consistent results, but the comparison will now be more 
demanding, given the smaller errors arising from current 
datasets. 



7 SUMMARY 

We have placed new constraints on the values of the basic 
cosmological parameters, using an up-to-date compilation of 
CMB data and the galaxy power spectrum measured from 
the final 2dFGRS by Cole et al. (2005). We have carried out 
a comprehensive exploration of parameter space, consider- 
ing five, six and seven parameter models, making different 
assumptions about the priors used for certain parameters. 
Our main results can be summarized as follows: 

• A model in which five parameters are allowed to vary 
does a remarkably good job of describing the currently avail- 
able CMB and LSS data. 

• There is an impressive level of agreement between the 
results obtained for CMB data alone and for CMB data plus 
the 2dFGRS power spectrum data. If the 2dFGRS P{k) 
is replaced by the SDSS P{k), there is some tension be- 
tween the parameter values preferred by the CMB and SDSS 
datasets. 

• For some parameters, for example the physical density 
of dark matter, Hubble's constant and the amplitude of den- 
sity fluctuations, there is a significant tightening of the al- 
lowed range of parameter space when the 2dFGRS P{k) is 
included in the analysis. In particular, we infer a density 
significantly below nm = 0.3. 

• We find some evidence for a departure from a scale 
invariant primordial spectrum of scalar fiuctuations. Our 
results for the scalar perturbation spectral index are only 
marginally consistent with the scale invariant value ns = 1; 



Cosmological parameters from CMB and 2dFGRS 21 



this spectrum is formally excluded at the 95% confidence 
level. However, this conclusion is weakened if we drop the 
assumption that the universe is fiat or allow neutrinos to 
have a mass. 

• We place new limits on the mass fraction of massive 
neutrinos: fu < 0.105 and ^ < 1.2 eV at the 95% level. 

• Several parameters are sensitive to the choice of prior 
for the optical depth to the last scattering surface, r. 

• We find that a wide range of closed universes are con- 
sistent with the CMB data. This range is restricted if we 
also consider the 2dFGRS P{k) data. If we further assume 
a prior of r < 0.3, then the preferred spatial curvature is 
close to fiat. 

• We confirm the evidence previously reported by Efs- 
tathiou et al. (2002) for a non-zero dark energy contribution 
to the energy-density of the universe. 

• We find a redshift-independent equation of state for the 
dark energy of tode = — 0.851q i^^, consistent with a cosmo- 
logical constant. 

• Inflationary models with a scalar field potential with a 
V{(f>) oc <f)* term are ruled out by our analysis. 

The final message of this analysis is that meaning- 
ful comparison of the parameter constraints from different 
studies requires a clear listing of the free parameters and 
their prior distributions. Although current datasets measure 
many parameter combinations extremely well, important de- 
generacies remain. As we have discussed, there is a good 
chance that current measurements may be poised on the 
brink of rejecting the simplest 5-parameter model in favour 
of something more complicated. However, even if this step 
is taken, it will require much work before any such deviation 
from the standard model could have a unique interpretation. 
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